Astronomy <fe Astrophysics manuscript no. 7204 


© ESO 2008 


September 16, 2008 





VLBI observations of Jupiter with the Initial Test Station of 
LOFAR and the Nancay Decametric Array 

A. Nigl^, P. Zarka^, J. Kuijpers^, H. Falcke'^, L. Bahren'^, and L. Denis^ 



OO 

o 
o 

(N 
D 



6 



> 

o 

(N 

G\ 
O 
OO 

O 



X 



^ Department of Astrophysics (IMAPP), Radboud University Nijmegen, 6535 ED Nijmegen, The Netherlands 
^ LESIA, CNRS/Observatoire de Paris, 92195 Meudon, France 
^ ASTRON, 7990 AA Dwingeloo, The Netherlands 

Received 31 January 2007 / Accepted 31 May 2007 

ABSTRACT 

Aims. To demonstrate and test the capability of the next generation of low-frequency radio telescopes to perform high 
resolution observations across intra-continental baselines. Jupiter's strong burst emission is used to perform broadband 
full signal cross-correlations on time intervals of up to hundreds of milliseconds. 

Methods. Broadband VLBI observations at about 20 MHz on a basehne of ^ 50000 wavelengths were performed 
to achieve arcsecond angular resolution, lofar's Initial Test Station (lofar/its. The Netherlands) and the Nangay 
Decametric Array (nda, France) digitize the measured electric field with 12 bit and 14 bit in a 40 MHz baseband. The 
fine structure in Jupiter's signal was used for data synchronization prior to correlation on the time-series data. 
Results. Strong emission from Jupiter was detected during snapshots of a few seconds and detailed features down to 
microsecond time-scales were identified in dynamic spectra. Correlations of Jupiter's burst emission returned strong 
fringes on 1 ms time-scales over channels as narrow as a hundred kilohertz bandwidth. 

Conclusions. Long baseline interferometry is confirmed at low frequencies, in spite of phase shifts introduced by variations 
in ionospheric propagation characteristics. Phase coherence was preserved over tens to hundreds of milliseconds with a 
baseline of ~ 700 km. No significant variation with time was found in the correlations and an estimate for the fringe 
visibility of 1, suggested that the source was not resolved. The upper limit on the source region size of Jupiter lo-B 
S-bursts corresponds to an angular resolution of ~ 3 arcseconds. Adding remote stations to the lofar network at 
baselines up to thousand kilometers will provide 10 times higher resolution down to an arcsecond. 

Key words. Instrumentation: detectors, high angular resolution. Methods: data analysis, observational. Techniques: high 
angular resolution, interferometric. Planets and satellites: Jupiter and lo 



1. Introduction 

A new wave of low-frequency telescopes are about to en- 
ter astronomy. LOFAR (low frequency array) is a new- 
generation giant digital telescope and already its proto- 
type (lofar/its) is capable of performing broadband and 
phase-sensitive interferometry with high angular, time and 
frequency resolution. In this paper, the nature of Jupiter's 
burst emission is studied at low frequencies (< 40 MHz) on 
a baseline of ~50000 wavelengths and an angular resolu- 
tion of 3 arcseconds. We want to test over which band- 
width and time-scales coherence in Jupiter burst signal is 
preserved or can be reconstructed, in spite of ionospheric 
propagation effects causing phase variation, Faraday rota- 
tion, dispersion and refraction. We want to demonstrate 
that high resolution observations can be performed by cross 
correlations of signals from remote stations being added to 
the LOFAR telescope at baselines up to a thousand kilome- 
ters. 

Long baseline interferometry of Jupiter's decametric 
emission (dam) has been performed since 1965 on base- 
lines up to more than 7000 km, setting an upper l imit on the 
source size of 0.05 arcseconds (jLvnch et al.lll972h . Research 



was mainly done in the sixties and seventies using narrow- 



tization (ISlee & Hieeinslll965llDulk et al. 


1967llBrown et al. 


19681 Dulkl 


19701 IStannard et al.l Il970l 


Lvnch et al.l 1X972 


Lynch et al. 


1976D. 



Send offprint requests to: A. Nigl 
Correspondence to: anigl@astro.ru.nl 



Jupiter's DAM is produced below 40 MHz, limited by 
the maximum electron gyrofrequency at Jupiter's surface. 
The lower limit for ground-based observations is set by the 
plasma frequency of the Earth ionosphere between 5 MHz 
and 10 MHz. The emission mechanism of Jupiter's dam, 
mod ulated by lo, is b elieved to be cyclotron maser emis- 
sion (IWu fc Led[l979h . The emission can be very intense 
up to 10^~^ Jy is beamed into a thin hollow cone with its 
axis parallel to the local magnetic field. The hollow cone 
has a half-angle of 60° to 90° relative to the local magnetic 
field direction and has a cone wall thickness of a few degrees 
(iKaiser et al.ll2000l ) . This emission can be detected on Earth 
when the line of sight of the observer coincides with the wall 
of the hollow cone. The emission arrives in noise storms 
called long bursts (L-bursts) and i n narrow milli second 
burst s called short bursts f S-bursts). iPhillips et al.l (1986, 
|1987[ ) found a first indication for source growth during L- 
bursts from smaller than 10 arcseconds up to 60 arcseconds. 
For S-bursts a complex microstructurc at typically ^--^40 s 
period and ~30 kHz bandwidth has been observed. These 
sub-burst wave packets have been attributed to clouds of 
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electr ons moving along the lo flux tube (ift) (jCarr fc RevesI 
fT999|). 



Dulkl (|1970[) put an upper limit of 400 km to 4000 km 



at Jupiter on the size of the DAM source at 34 MHz, which 
corresponds to an angular resolution of 0.1 to 1 . S-burst 
analysis provides an upper limi t of about 20 km on the 
instantaneously emitting source (|Zarkalll996l) . which trans- 
lates to 6 milli-arcseconds at the distance of 4.25 AU. 

Jupiter's DAM emission is strongly elliptically polarized. 
The signal at frequencies near 20 MH z is on average 49% 
circularly and 87% linearly p olarized (jShaposhnikov et all 
ll997llQueinnec fc Zarkal l2001'). 

For the following analysis, Jupiter was observed with 
the two radio telescopes lofar/its (Initial Test Station 
of LOFAR, The Netherlands) and nda (Nangay Decametric 
Array, France), which arc described in Sect. ^ In Sect. 
21 correlations of simultaneous data recorded with ITS 
and NDA are described and then discussed in Sect. [S] 
Correlations were performed on dynamic spectra and on 
the time domain data. 



2. Radio telescopes 

2.1. ITS 

lofar's Initial Test Station was located at the core site of 
LOFAR near the village of Exloo, The Netherlands and con- 
sists of 60 inverted- V-shaped dipole antennae distributed 
over five spiral arms, its antennae work in the lowest lofar 
frequency band from 5 MHz to 35 MHz, which matches the 
emission band of Jupitcr-Io emission. 

The final lofar telescope will eventually consist of 
~15000 antennae, split up between a dense core of 32 sta- 
tions and 45 remote stations with 192 antennae each. These 
stations will be spaced on exponentially expanding shells. 
The antennae within the station will be placed on expo- 
nential shells with a maximum diameter of 60 meters. The 
total bandwidth from 30 MHz to 240 MHz wiU be covered 
by two different antenna types: a low band antenna (lba) 
for frequencies below 80 MHz, and a high band antenna 
(hba) from 110 MHz to 240 MHz. Below 30 MHz the gain 
of the low band antennae is reduced, but a fixed hardware 
limit is at 10 MHz. The lofar electronics is designed to 
accommodate a third type of antenna, which could cover 
frequencies below 30 MHz to have a broader bandwidth for 
Jupiter observations. With a maximum baseline of ~100 
km between stations, all these antennae will provide lofar 
with a 2 to 20 arcseconds angular resolution and millijan- 
sky sensitivity, dependent on frequency. Each dipole's signal 
will be amplified in an l na (low-noise amp lifier) , filtered in 
frequency and digitized (jFalcke et al.ll2006l ). More detail on 
LOFAR and its test stations can be found at|www . lofar . org[ 
Additionally, a formulation of a science case with long base- 
lines by ext remely remote LOFAR stations is described by 
rVo^(l200l . 

In the case of its, high performance analog-to-digital 
converters (adcs) with a dynamic range of 12 bits, supply 
a rate of 80 million samples per second. The sample rate 
results in a time resolution of 12.5 ns and a data rate of 
~160 MB per second and dipole. Four of those ADCs are 
mounted on two TiM-boards (twin-input module) which arc 
accommodated in each of the 16 standard PCs in the data 
acquisition (daq) container. The TiM-boards can handle 
two digital data streams and store each stream in 1 GB of 



RAM. This amount of RAM limits the DAQ to 6.7 s. From 
2005 on, 30 antennae measured dual polarization, since two 
orthogonal dipoles were mounted in each antenna structure. 
Each antenna detects the whole sky with a primary beam 
of about 90° full width half maximum (fwhm) centered on 
the zenith. Using digital beam-forming the array can be 
simultaneously pointed at several sources in the sky. For 
its the digital beam-forming was performed off-line and is 
described in Sect. [3l The absolute timing of its was up- 
dated with Network Time Protocol (ntp). The beam size 
of ITS at the zenith results in an angular resolution of 3, at 
a wavelength of 10 meters (30 MHz) and for a maximum 
baseline of 180 meters. 

2.2. NDA 



www . obs-nancay . f r/dami 

4 antennae. Each antenna 



The Nangay Decametric Array 
is a phased array made up of 1 
consists of eight conducting wires (monopoles) wound up 
on vertical support cables on a conical surface. Each op- 
posite pair of monopoles is connected as a spiral dipole. 
The array is split in two groups of 72 left-handed and 72 
right-handed wound antennae to measure circularly polar- 
ized radio waves. All 72 antennae of one cluster measure 
the same polarization direction and can be phased up with 
an analog b eam- former. A f ull de scription of the array can 
be found in lBoischot et all (|1980[ ). 

The receiver used for the observations in 2005 provided 
the same type of data as its by sampling in the first Nyquist 
zone with 80 MHz over a frequency band from 8 MHz to 
40 MHz. A DAQ board (Signatec PDA14) was used for the 
digitization with a 14-bit ADC inserted in a powerful PC, 
which allows continuous DAQ of a stream of 80 million sam- 
ples per second. These data are written as short integers 
(16-bit) to disc. The receiver allows full signal detection 
below 40 MHz (a low-pass filter at 40 MHz and high-pass 
filter at 8 MHz, 16 MHz or 24 MHz are used) without use of 
a local oscillator for frequency down conversion. Absolute 
timing of the instrument is provided by a receiver of the 
Global Positioning System (gps). 

3. Observations 

Triggered by NDA, we took simultaneous datasets dur- 
ing periods of Jupiter burst emission. Simultaneous 
Jupiter snapshots of 6.7 s (1 GB) at its and 20 s 
(~3 GB) at NDA were recorded. The trigger was man- 
ually executed at NDA while monitoring the real-time 
display of NDA for actual burst activity. The real- 
time monitor produces and displays dynamic spectra 
at a rate of one spectrum per second. The monitor 
can also be found online with a 30 second refreshment 



time ( www . obs-nancay . f r/dain/daiii_visu . html) . The day 
of observation was selected from NDA's monthly prob- 
ability plots ([www ■obs-nancay.fr/dain/daiii_proba .htm 
iLeblanc et allliggsf r 

Both telescopes recorded the full signal during overlap- 
ping periods. The data taken with NDA are available in dig- 
ital format and contain the digitized output of the analog- 
formed beam in the direction of Jupiter. For its the raw 
digitized data for all antennae were recorded. Total inten- 
sity was not calculated due to uncertain relative gain cali- 
bration. Since a single dipole antenna measured the whole 
sky at the same time, it detected also the sky background 
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as well as strong radio frequency interference (rfi). The 
resulting high noise level makes it impossible to observe 
Jupiter's emission in one antenna with enough contrast. 
The phasing of the antennae by digital beam-forming, in the 
direction of Jupiter, was performed to increase the signal- 
to-noise ratio by suppressing the RFI coming from the hori- 
zon. 

To beam-form and display the data, the following steps 
were performed. First, a Fourier transform, however the 
computation time of the fast Fourier transform (fft) in- 
creases as N logN, where N is the number of samples. 
Thus, splitting the antenna data in time-blocks is more ef- 
ficient than processing one long time-series. However, the 
splitting requires extra bookkeeping when reconstructing 
the beam-formed time-series needed for displaying and for 
signal correlation. 

Second, a Hamming window was applied to the raw data 
blocks to keep the leakage in frequency domain into neigh- 
bouring channels as low as possible. Third, the fft was 
applied on blocks of 16384 samples. Fourth, the interesting 
band or bands were selected from the frequency domain. 
Fifth, the beam-forming phases were applied to each fre- 
quency channel. The beam-forming phases correct the dif- 
ferential light travel time to each dipole-element relative 
to a reference position or antenna. They are calculated for 
the exact fractional sample delay, that shifts the antenna 
data for the exact desired beam direction. The chosen block 
size determined a time resolution of 0.2 ms and a spectral 
resolution of 5 kHz. After averaging all the phase-shifted 
antenna spectra the resulting beam-spectra were stacked 
for each block and displayed in a dynamic spectrum. The 
process of true time-delay beam steering by phase-shifting 
and element averaging on ADC data is called digital beam- 
forming. 

For the reconstruction of the time-series, the beam- 
formed blocks had to partly overlap. The start and end 
of each beam-formed block will not contain data from all 
antennae, due to the delays applied to each antenna during 
beam-forming. For the resulting beam-formed time-series, 
ultimately only intervals over which delay-shifted data from 
all antennae exist are usable. 

For visual comparison of the simultaneous data, the dy- 
namic power spectra from the two telescopes are plotted in 
Fig. [TJ The 48 antennae of its and the 72 antennae of NDA 
were beam-formed in the direction of Jupiter. The spec- 
tra were computed with 0.2 ms time resolution and conse- 
quently a spectral resolution of 4.9 kHz. The Jupiter sig- 
nal received was expected between 20 MHz and 30 MHz, 
where only a few narrow band lines contaminate the spec- 
tra. Below 20 MHz strong RFI is dominating the spectra 
and above 30 MHz the bandpass filter of its rolls off. The 
Jupiter-Io S-burst emission is clearly visible between 20.7 
MHz to 23.9 MHz and several identical features can be iden- 
tified. The emission drifts in frequency, which translates to 
a movement of the source region along Jupiter's magnetic 
field lines. These S-bursts last for a fraction of a second 
down to a few milliseconds. The horizontal bands in the 
ITS plot are due to Faraday rotation. The power spectra of 
the two perpendicular polarization directions of its show 
strong Faraday fringes (see Fig. [2]). The spectra have a spec- 
tral resolution of 1.2 kHz and are integrated over 6.7 s. The 
plotted frequency range from 21.7 MHz to 22.95 MHz covers 
the strong Jupiter burst emission and contains only little 
RFI (narrow vertical lines). The Faraday fringes were fitted 
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Fig. 1. Dynamic spectra of a simultaneous observation with 
the Nangay Decametric Array (nda) [top] and lofar's 
Initial Test Station (lofar/its) [center], on November 
30, 2005 at 07:15:12 UTC for a duration of --3.5 s. Jovian 
S-bursts show up as the discrete structures mainly con- 
fined to the band 20.7 MHz to 23.9 MHz. Dark horizontal 
Faraday fringes are superimposed on emission observed at 
ITS. Vertical lines in the NDA spectrum are 5 ms pulses in- 
tentionally added every second to the data for timing pur- 
poses. Horizontal lines are caused by radio frequency inter- 
ference (rfi). The third panel [bottom] shows the result- 
ing cross correlation coefficients versus time and frequency, 
computed from data blocks of 2.5 ms duration and ^--150 
kHz bandwidth. 



with a sinusoidal function according to Faraday's law. The 
resulting rotation measure of RM = (1.15 ± 0.07) rad m~^ 
is in good agreement with the rotation measure determined 
from total electron content values provided by the GPS. 
For the East-West polarization direction, the maxima of 
the Faraday fringes are indicated with dots. The width of 
the Faraday fringes was measured as 100 kHz to 150 kHz. 
Around the frequencies of the eight maxima, the signal was 
filtered for the cross correlations described in Sect. 14.31 For 
a more detaile d analysis of the F araday rotation and previ- 
ous results see iNigl et aD (|2005[ ). 

4. Long baseline interferometry 

For this work, Jupiter was observed with its and NDA on 
a baseline D of 702 km, corresponding to an angular res- 
olution 6 of 3.3 at a frequency of = 30 MHz. For this 
resolution, Jupiter burst emission regions can be consid- 
ered as point sources. The smallest structure that could 
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Fig. 2. Spectra of the two linear polarizations of lofar's 
Initial Test Station (lofar/its) [red: E-W antennae dipolc 
beam and blue: N-S antennae dipole beam]. The best-fit 
sine function (in frequency, not intensity) of the computed 
Faraday angle is superimposed to both ITS spectra. For 
the East- West polarization direction the maxima of the 
Faraday fringes are indicated with dots. 



be resolved on Jupiter at the time of the observation was 
13600 km in diameter, which is four times the size of lo. 



4.1. 2D Cross correlation of dynamic spectra 

The synchronization of the absolute timing of its and NDA 
data was of the order of one second. In order to restrict the 
size of the time window for correlation, a normalized 2D 
cross correlation of the dynamic spectra shown in Fig. [1] was 
performed. This cross correlation provided a coarse fringe 
search for the time delay and an offset in frequency between 
the two datasets, which is called a clock model in VLBI ter- 
minology. Before correlation, the dynamic spectra have to 
be generated with high time resolution for measuring the 
time delay, and with high spectral resolution for determin- 
ing the frequency offset. The higher the resolution in time, 
the lower in frequency, maintaining a time-bandwidth prod- 
uct of unity {dv x dt = 1). Thus for high time resolution, 
fine spectral structure is lost, which reduces the quality of 
the subsequent 2D cross correlation, but gives access to an 
accurate time delay. An offset in frequency can be found by 
choosing high spectral resolution. 

For the observation performed on November 30, 2005 at 
07:15:12 UTC, the result of a 2D cross correlation, of two 
slices of the dynamic spectra from its and NDA is given in 
Fig. [21 In this correlation window, maximum correlation is 
found at zero frequency shift and at a time delay of 4.9 ms. 
The clock offset resulted in 10.0723328 s with a resolution 
of 1.6 fis (its started to record later than nda). This clock 
offset (clock model) serves as input for the time-scries cross 
correlation. Additionally, a shift in frequency of 191 Hz with 
a resolution of 38 Hz was obtained, which will be explained 
in Sect. l4.5l as being mainly due to a difference between the 
sample rate of both telescopes, caused by the deviation of 
the ITS system from the nominal sample rate. 

The sample rate was measured for NDA with (79999998± 
0.5) samples per second differing only slightly from the 
nominal sample rate of 80 megasamples per second. This 
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Fig. 3. Correlation coefficient matrix for two dynamic 
spectra slices. The intensity in the image displays the cross 
correlation coefficient labeled with the color wedge on the 
right side of the plot. 



difference causes an offset from the expected and recon- 
structed frequency in the dynamic spectrum of -1-0.6 Hz (for 
NDA at 22.3 MHz). Furthermore, due to Earth rotation, the 
radial speed of the two telescopes with respect to Jupiter at 
the time of the observation differed by 38 m s~^, with its 
moving towards Jupiter as seen from NDA. Therefore the its 
signal was stronger Doppler shifted to higher frequencies by 
2.8 Hz. Both frequency offsets lie within the frequency offset 
resulting from the 2D cross correlation above. 



4.2. Coherence time and coherence bandwidth 

Since emission from Jupiter arrives in bursts, the best 
signal-to-nose ratio is achieved by correlating with integra- 
tion time and frequency resolution matched to the dura- 
tion and bandwidth of the bursts. The emitting electron 
bunches at Jupiter are traveling along the ift at several 
thousands to several tens of thousands of kilometers per sec- 
ond. The emission is received as subpulses (wave packets) 
with about 40 fj,s to 90 fiS duration and a m inimum band- 
width of 30 kHz. as it has been observed bv lCarr fc Reved 
(1999,). Therefore the blocks for correlation should be much 
longer than a wave packet of 90 /iS (thus ~ 1 ms), and have 
a frequency resolution coarser than 30 kHz, but not more 
than the Faraday fringe spacing of 150 kHz. The upper limit 
for the correlation time window is dependent on the time- 
scales of changes in the Earth ionosphere and has yet to 
be identified. The discussed limits constrain the optimum 
window for the cross correlation. 
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Time Lag [ms] 

Fig. 4. Result of the time-series signal cross correlation of 
the NDA and its data sets of the Jupiter observation on 
30.11.2005 at 07:15:12 UTC. The (red) peaked curve in the 
plot was obtained by averaging about 1000 correlated 1 ms 
time-series blocks containing Jupiter burst signal at eight 
bands of 75 kHz bandwidth. The (blue) flat curve was ob- 
tained on the same number of blocks without emission. The 
peak maximum lies 210 sigma above the noise level with 1 
s effective integration. 



4.3. Cross correlation of time series 

The final cross correlation was done in the frequency do- 
main following Parseval's Theorem. The fringe amplitude 
was calculated with the so-called fringe correlation method 
as the square root of the su m of the squares o f the quadra- 
ture correlations (compare iLvnch et al.lll976l ). The cosine 
component was obtained by a second cross correlation with 
a quarter-period shift at the center frequency of the se- 
lected subband introduced to one of the time series. The 
quarter-period shift was applied in the frequency domain 
as a phase-gradient, which enables shifts by a fraction of a 
time sample. 

We used the cross correlation theorem, which performs 
the cross correlation of one block, for all possible lags, by 
only three Fourier transformations and one multiplication 
of the two data arrays. The lags in the cross correlation 
result in sample steps of 12.5 ns. 

Since the observed Jupiter-Io emission is not continu- 
ous, noise gating was applied and only time blocks with 
received power exceeding the noise-fioor were added after 
cross correlation according to the fringe correlation method. 
Integrating about one thousand 1 ms blocks containing 
burst emission at eight bands of 75 kHz bandwidth, re- 
sulted in an average cross correlation peak of 210 sigma 
significance compared to the same number of integrated 
blocks containing only background noise (see Fig. 2]). 

In total, narrow bands with bandwidths of (30, 75 and 
150) kHz at center frequencies of (21.73, 21.88, 22.04, 22.20, 
22.36, 22.53, 22.70 and 22.88) MHz were analyzed. These 
frequencies were chosen as they lie at the strong maxima of 
the Faraday fringes from the its power spectra, which are 
shown in Fig. [21 

The bottom panel in Fig.[l]displays the cross-correlation 
coefficients in a dynamic spectrum for 2 seconds between 



20.6 and 23.8 MHz, computed from 2.5 milliseconds time 
blocks over all bright Faraday fringes. 

The distribution of cross correlation coefficients of 500 
correlated 1 millisecond time blocks (selected by noise gat- 
ing), was maximum for a bandwidth of 30 kHz at a center 
frequency of 22.53 MHz resulting in C = 0.6 ± 0.1 (one 
sigma statistical error stated). 

4.4. Fringe visibilities 

The cross correlation coefficients of the time blocks are 
proportional to the fringe visibilities, which ultimately 
car ry information on th e source size, following for exam- 
ple |PhiiiiiiIiES ([l98l : 



7s 



Gr 



kCrji 



' [Psi + Pbi){Ps2 + Pb2) 



PS1PS2 



(1) 



Here 7s is the instrument-corrected fringe visibility due 
to the source structure; Gr is the fringe washing factor at 
the residual geometrical delay r; C is the cross correlation 
coefficient; k and 77 are instrumental decorrelation coeffi- 
cients; Ps is the power from the source and Pb the power 
from the Galactic background both recorded at the two 
telescopes. The instrumental decorrelation coefficient rj ac- 
counts for any phase behavior of filters, amplifiers, cables 
and connectors of both telescopes together. For its/nda 
this factor was assumed to be maximum 15% and thus 
7? = 1.15. 

Since only circular polarized parts of the Jupiter sig- 
nal were measured by NDA and correlated with only linear 
polarized components of its, a maximum cross correlation 
coefficient of Cmax ~ \/2/2 can be achieved. Therefore the 
second instrumental coefficient was adopted to fc = ^/2. 

The fringe washing functio n, Gt, is the Fourier t rans- 
form of the receiver passband ([Thompson '^MMbEB)- For 
a rectangular passband, Gr is a sine function: Gt{t) = 
sinc{T: Avt) . Here, Av is the bandwidth of the correlated 
signal and r is a residual geometrical delay, which is the 
difference between the true delay and the actually applied 
delay for the correlation: t ~ Tg ~ ti. In the off-line corre- 
lations performed for this paper, the lag of the maximum 
correlation was determined for each block, which resulted 
in T = and G = 1. The two terms [signal -\- noise) / signal 
under the square root in Eq. [T]are estimated to be ~ 1.04 
for NDA and ^ 1.06 for ITS. A global background correc- 
tion factor of 1.1 is obtained. The resulting visibilities 75 
are consistent with one, and thus with unresolved S-burst 
sources, as expected. 



4.5. Cross correlation lags 

The lag in which the correlation maxima laid, and hence 
the signal arrival time difference between the two datasets, 
changed by (—8.8 ± 0.1) fjLS (one sigma statistical er- 
ror stated). The minus sign means that the signal arrived 
progressively earlier at ITS. This drift in the lag is mainly 
caused by a constant relative drift of the clock quartz of 
ITS with an accuracy of AF/F ~ 10^^. Here F is the 
nominal frequency and AF the maximum deviation in- 
cluding errors due to calibration, temperature variation, 
shock, vibration and aging over 10 years. The accuracy 
of the NDA clock is much better and was measured by 
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GPS with AF/F ~ 1.25 x 10~^ and thus accounts for a 
maximum drift rate of only 12.5 ns s~^. The motion of 
the source with an average paraUcl electron velocity of 
- 22000 km s^^ (|Zarkalll996l ) translates to a drift of only 
—367.82 ns s"-^ at the Earth. Interplanetary multipath scin- 
tillation can cause delays of the order of only a few nanosec - 
onds assuming a thin scattering sheet (jPhillips et al.|[l987[ ). 
The Earth rotation during the short observation only ac- 
counts for -1-125.17 ns s~^ and ionospheric effects change 
on time-scales longer than the duration of our observation 
and therefore were approximated by a differential change 
in the local slant tec values, resulting in a contribution of 
-0.032 ns s-i. 

Additionally, the drift of the its clock causes a deviation 
between the sample rates of NDA and its of 704 Hz. The big- 
ger the time- window for the cross correlation, the larger the 
deviation in time per correlated sample gets, which reduces 
the quality of the correlation. Furthermore, this difference 
in the sample rates translates into an offset in frequency 
of 196 Hz between the dynamic spectra of both telescopes, 
which is consistent with the frequency offset of 191 Hz ob- 
tained from the 2D cross correlation of the dynamic spectra 
with a resolution of 38 Hz (Sect. 14. ip . 

5. Conclusions & outlook 

Using the prototype its of the new generation radio tele- 
scope LOFAR and the Nangay Decametric Array, strong 
lo-induced Jupiter radio emission, including S-bursts, was 
detected over 6.7 seconds with only dozens of antennae. 
The beamed dynamic spectra obtained with its and NDA 
show excellent agreement (see Fig. [T]) . High time and fre- 
quency resolution allowed us to determine offsets by two- 
dimensional cross correlation of the dynamic spectra with 
a precision of microsecond accuracy in time and hertz ac- 
curacy in frequency. The latter method revealed a shift in 
frequency, due to the relative radial speed and an offset in 
sample frequency between the two telescopes. 

The baseline of the interferometer of 702 km provided an 
angular resolution of three arcseconds. Since the emitting 
source region is believed to be smaller than 0.02 arcseconds 
(|Zarkalll996f ). the source region of the received emission re- 
mains unresolved, as confirmed by the obtained visibilities 
being close to one. The identified delays at the maxima 
of the time-series cross correlations showed a linear drift, 
caused by a large clock rate at its. Correcting for this con- 
stant drift, no significant variation was found within the ob- 
servation of 6.7 s. One second effective integration in time 
and several hundred kilohertz bands in frequency resulted 
in an 210 sigma detection above the noise level. 

This is the first time that VLBi-type observations have 
been performed over a broad relative bandwidth with base- 
band digitization, and direct correlation of signals with a 
large dynamic range. Cross correlation fringes on the beam- 
formed time-series were found on time-scales of millisec- 
onds. As demonstrated with the its prototype, lofar will 
probably reach arcsccond angular resolution, lofar will be 
pioneering in being capable of performing VLBi-type obser- 
vations as described in this article with millijansky sensi- 
tivity. The present results support the European extension 
of LOFAR with very remote stations up to several hundred 
kilometers. A baseline of a thousand kilometers will serve 
LOFAR with a spatial resolution of an arcsecond at 30 MHz. 
Such a resolution will enable detailed studies of the Jupiter 



magnetosphere as described in IZarkal (f2004l ) among many 
other astrophysical applications ( Vogtll2006[ ). 

A follow up of these long-baseline interferomctry studies 
is planned with NDA and lofar/csI (Core Station 1). The 
follow up experiments will be performed to quantify at what 
duty cycle time-phase-coherence is preserved, as a function 
of time and frequency. Furthermore, weaker radio sources 
may be used in addition to Jupiter. 
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